Dynamics and kinetics of quasiparticle decay in a nearly-one-dimensional degenerate 

Bose gas 



O 

(N 

> 
O 

in 



<Z3 

-t— > 



i 

C 

O 
o 



> 

m 

ON 

m 
O 



X 



I. E. Mazets 

Vienna Center for Quantum Science and Technology, Atommstitut, TU Wien, 1020 Vienna, Austria; 
loffe Physico-Technical Institute, 194021 St. Petersburg, Russia 

(Dated: ) 

We consider decay of a quasiparticle in a nearly-one-dimensional quasicondensate of trapped 
atoms, where virtual excitations of transverse modes break down one-dimensionality and integra- 
bility, giving rise to effective three-body elastic collisions. We calculate the matrix element for the 
process that involves one incoming quasiparticle and three outgoing quasiparticles. Scattering that 
involves low-frequency modes with high thermal population results in a diffusive dynamics of a 
bunch of quasiparticles created in the system. 

PACS numbers: 03.75.Kk,05.30. Jp 



I. INTRODUCTION 

An uniform one-dimensional (ID) system of indistin- 
guishable bosons interacting with each other via pairwise 
delta-functional potential is known to be integrable and 
described by the Lieb-Liniger model [1]. In an integrable 
system, the number of integrals of motion equals to the 
number of degrees of freedom. Such an equality, on the 
one hand, facilitates analytical treatment of integrable 
models and, on the other hand, make their dynamics 
quite distinct from the dynamics of non- integrable many- 
body systems. In the course of its evolution, an inte- 
grable system always "remembers" the information on 
its initial conditions, and relaxation to a thermal equilib- 
rium state does not occur. A general review of integrable 
models can be found, e.g., in Ref. Q. The Lieb-Liniger 
model can be implemented in an ultracold-atom exper- 
iment in optical lattices Q or on atom chips [4]. The 
conditions of the ID regime require both the tempera- 
ture and mean interaction energy per atom being well 
below the excitation quantum huj± of the radial trapping 
(harmonic oscillator) Hamiltonian: 



T < hcoi 



5lD«lD iS ?1W_L- 



(1) 



Here temperature T is measured in units of energy (i.e., 
we set &b = 1) and mo is the linear density of atoms. 
Since the effective ID interaction strength for atoms 
under radial harmonic confinement is (far below the 
confinement-induced resonance Q) gxu = 2hui±a s , where 
a s is the atomic s-wave scattering length Under these 
conditions [6j, the latter of two Eqs. (JTJ reads 



ni D a. 



< 1. 



(2) 



If Eqs. (p} hold then radial motion of atoms is confined 
to the ground state of the radial trapping Hamiltonian. 
Ultracold atomic systems in the ID regime can be pre- 
pared in optical lattices [7| and on atom chips [8|- How- 
ever, in reality no system is perfectly ID, but the actual 
question is, on which timescale it can be described as 
ID. Since the thermal population of the excited states 
is strongly suppressed by the exponential Boltzmannian 



factor, already T s» 0.2 huj± [8| corresponds to a deeply 
ID regime in terms of thermal excitations. On the other 
hand, the influence of atomic interactions is a more com- 
plicated and interesting issue. For ultracold atoms under 
tight lateral confinement one-dimensionality and, hence, 
integrability are lifted by atomic interactions causing vir- 
tual population of excited radial modes. The role of the 
virtual radial excitations in the dynamics of ultracold 
atomic gases in tight waveguides has been first studied 
in the context of macroscopic flow of de gen erate atomic 
gas through a waveguide [9| and decay [l(| or inelastic 
collisions of mean-field solitons. On the microscopic 
level, a second-order perturbative collisional process with 
a radially excited virtual (intermediate) state give rises to 
effective three-body elastic scattering that has been sug- 
gested as the source of thermalization in ultracold atomic 
gases on atom chips [H, EH • 
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FIG. 1: Feynman diagrams of the processes studied in the 
present paper: (a) splitting of an elementary excitation into 
three excitations, (b) the reciprocal process. 

Damping of a fast particle motion in a quasi-lD 
bosonic system due to the effective three-body process 
shown in Fig. [1] has been studied recently [Hj], and the 
contribution of events with small momentum transfer to 
the damping rate has been found significant in a certain 
parameter range. The reason is the bosonic amplification 
of scattering into modes with high thermal population 
14] . In the present paper we derive (i) the vertex for the 
diagrams in Fig. [1] and (ii) a kinetic equation describ- 
ing damping relaxation of a quasiparticle in a degenerate 
nearly-lD bosonic system and its diffusion-type limit de- 
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scribing the elementary excittation dynamics induced by 
collisions with a small momentum transfer. 

Using our theoretical estimations, we conclude that 
the three-body elastic processes responsible for the in- 
tegrability breakdown can be detected experimentally by 
observing the damping dynamics of an ensemble of fast 
particles on relatively short time scales. 



where K 3b = \£,n m a s . 

We expand the phase and density operators in plane 
waves. The average ID density tt-id = N/L is an integral 
of motion has a well-defined value since the total number 
N of atoms is fixed. The expansions read as 



rim 



= *n 



k 



II. PHYSICAL MODEL 

The Hamiltonian of a ID system of identical bosons 
with an additional term accounting for effective three- 
body interactions via virtual radial excitations reads 



H = H LL + H 3 b, 



(3) 



where is the Hamiltonian of the Lieb-Liniger model 
written in terms of the ID number density fi and phase 
4> operators [H, Il6j 



#LL = 



dz 
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(4) 

where m is the atomic mass and L is the quantization 
length. The commutation rule for the phase and density 
operators is [(f>(z), n(z')] = —id(z — z'). The three-body 
collisional dynamics is taken into account by the term 



dz n 3 , 



(5) 



where £ = 4 In | . Note that an error in the prefac- 
tor of H ni [l2j has been corrected later [H, 0. Of 
course, Eq. (|5|) is only a lowest term in expansion of the 
integrability-breaking interaction in powers of the small 
parameter nrDa s . The full effective interaction due to ex- 
citation of radial modes contains cubic, quartic, quintic 
in h terms and supports a stable ground state. This can 
be seen from the analogy with the mean-field variational 

approach [lEl. However » if Ec l- © holds then Ec L- © 
is sufficient for perturbative calculation of the rates of 
the processes shown schematically in Fig. [TJ 

From now on, we scale energy to <7id^id and set the 
healing length (h — h/ ^/mgiD^-iD as the length unit, 

i.e., H = gmniDH and z — C^z. In what follows, we 
omit, for the sake of compactness of notation, the bar 
over scaled values. It will not lead to confusion, since we 
will mention explicitly the return to usual units. In the 
new, dimensionless form Eqs. (|4l [5]) read simply as 



The global phase $o, subject to phase diffusion even if 
-/V = const [17j . will not appear in further expressions. 
The sums in Eq. are taken of non-zero (both positive 
and negative) integer multiples of 2tt/L. The commu- 
tation rule for the phase and density operators in the 
momentum representation are 



(9) 



By re-grouping terms we express the full Hamiltonian 
of the system as H = E gI + + H a + H n i, where E gI 
is the ground state energy, the two next terms represent 
the harmonic 



H h = 



k 



n 1D k 2 



A' 2 



G T 
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(10) 



and anharmonic 
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(ii) 

parts of the Lieb-Liniger model Hamiltonian. Operator 
Vk in Eq. (fTTj) is defined as 



"ID 



dz e 
1 



sfh; 



-ikz 



E 



'ID 



h k _ k ,h k , + 



(12) 



G rn = 1 — |£ 71 iDa s in Eq. ([TU]) is the effective strength 
of pairwise interactions renormalized by the presence of 
the three-body collisions. In usual units, the coupling 
strength giD renormalizes to gJ D = Ggiu- However, since 
Eq. ([2j holds by assumption of the ID regime, we use in 
what follows an approximation G rn ~ 1. 

The three-body interaction term that makes our model 
non-integrable is 



K 3 b 



\fhn 



2 ^'"-1-1 
ID q ,q> 



(13) 



H LL = I dz 



1 d(j) „ d(j> \ „_ x ( dh 
2~d~z n ~d~z + t 1 \~d~z 
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Neglect for a moment the anharmonic part ([TT]) of the 
Lieb-Liniger model Hamiltonian. The retained harmonic 
part (fTOj) is easily diagonalized [15]: — J2k e kb k b k , 
where Bogoliubov elementary excitation energy is 



e k = |fc|Vl + fc 2 /4, 



(14) 
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and bk (b\) is the annihilation (creation, respectively) 
operator for an excitation with the wavenumber k. The 
density and phase operators are then 



where 



\k\ 



2e fc 



. (15) 



(16) 



is the static structure factor of the quasicondensate at 
zero temperature. In this (harmonic) approximation the 
operator fik annihilates one elementary excitation with 
the wave number k (or creates an excitation with the 
wave number —k), and the interaction (fTS"]) can cause 
only conventional Beliaev [H, HH or Landau [ill Gil US] 
damping. However, these types of relaxation caused by 
decay of an elementary excitation into two excitations 
or by the reciprocal process are completely suppressed in 
ID by energy and momentum conservation, in contrast to 
the 2D and 3D cases. To calculate the rates of processes 
shown in Fig. [TJ we need to diagonalize (in some approx- 
imation) the whole Lieb-Liniger Hamiltonian + H a . 
The unitary transformation providing this diagonaliza- 
tion transforms h q into an operator containing correction 
term that is nonlinear in creation (annihilation) opera- 
tors of elementary excitations and thus allows for the 
processes shown in in Fig. [T] The transformation is very 
similar to the polaronic transformation [2l|, but, unlike 
the latter, does not involve impurity particles. In other 
words, the idea is to demonstrate that an elementary 
excitation in the Lieb-Liniger model corresponds to a 
phase-density wave at a certain wavelength, dressed by 
virtual phase-density waves. 

To obtain analytic results, we set Vk ~ fik and thus 
obtain 



ft ~ 1 V w ' 



1,1 



i 



IJ q"-q-q' Vq' 



4n? D 



71 — q — q' TvqTlrql 



(17) 

Before dealing with the particular Hamiltonian of the 
problem, we discuss the diagonalization procedure in gen- 
eral. 



A. Overview of the diagonalization procedure 

Consider a Hamiltonian 

R=R +eW, (18) 

where the unperturbed Hamiltonian Rq is diagonal in the 
orthonormal basis {|A)}: 



H = ^4 0) |A)(A| 



(19) 



The perturbation operator contains the small parameter 



E 



(0) 



E 



(0) 



In what 



e. More precisely, e(A'|W|A) < 

follows we assume that the diagonal matrix elements of 
>V is the basis {|A)} are zero (if not, we can eliminate 

them by including e(A|VV|A) into ). 

The basis functions |A) that diagonalize the full Hamil- 
tonian jiHii 



H = 5> A |A)(A| 



(20) 



are related to the old basis via unitary transformation 

|A)=e*|A), (21) 

where the operator R = — BJ is anti-Hcrmitian. In the 
pertzurbative approach, R can be expanded in series in 
the powers of the small parameter 



R = ^Rj 
i=i 



(22) 



The standard quantum-mechanical perturbation theory 
[22] ] yields explicit expressions for the lowest-order terms 
in the expansion (|22p . in particular, 



|A')(A'|W|A)(A| 

E (o) _ E (o) ' 



(23) 



Instead of transforming wave functions, we can trans- 
form operators: 



A 



e- R Ae R 



= A-[R,A] + -[R,[R,A}] 



(24) 



Here A stands for an arbitrary operator. If, for example, 
we substitute Ti instead of A into Eq. ([24]) and apply 
Eq. we obtain 



R — Rq ~ £ 2 A"H2 



(25) 



where AR 2 = [R 2 ,R } + [Ri,W] - ±[Ri,[Ri,Hq]] is 
the term describing second-order correction to the un- 
perturbed energies and, in a case of coupling to con- 
tinuum, widths (to the second-order approximation) of 
decaying states. The term linear in e is absent because 
(A|W|A) = for all A by assumption. However, the sig- 
nificance of Eq. ([Ml) transcends far beyond derivation of 
eigenenergies E\. Eq. (|24p provides the transformation 
rule for any arbitrary operator, in particular, the atomic 
ID density operator. The first-order approximation for 
Eq. ([24]) reads as 



(26) 



A « A- e[Ri,A], 
where R± is given by Eq. 
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B. Application to the Lieb-Liniger model 
perturbed by three-body collisions in the 
weakly-interacting regime 

We take now Hq = Hh and eVV = H a , where the an- 
harmonic part of the Hamiltonian is approximated by Eq. 
(fTT)) . Direct calculation of the corresponding R as eR\ 
is rather involved, therefore we use the following ap- 
proach. First, we note that R\ is cubic in creation (anni- 
hilation) operators of elementary excitations, being the 
eigenmodes of H^. This means that the general form of 
eR\ is also cubic in terms of "bare" phase and density 
operators: 

E Ri= iii =^7fY. ( B ~ q ~ q ' q *M- q -A* + 

9,9' 

B q _ q q '_ q ,4> q n- q - q 4 q ' + B~l~ q ng4>-g- q >h q , + 

B—q — q' q q /7l q Tl—q—q'Tl q 'J. (27) 

The summation in Eq. (|27[) is taken over non-zero 
wavenumbers, i.e., if one of the wavenumbers q, q' or q+q' 
equals to 0, then the corresponding term is dropped off 
the sum. 

It is easy to prove that the general solution of an opera- 
tor equation [Q, Hq] = eVV, which follows from Eq. ([25]). 
is Q = eR\ + X, where I is any operator commuting with 
the unperturbed Hamiltonian, [I, Ho], i.e., correspond- 



ing to a certain integral of motion. However, because of 
energy and momentum conservation in ID, there is no 
operator, which has the form of Eq. ([2~T1) and commutes 
with _ffh- In other words, the sought for polaronic unitary 
transformation can be uniquely determined (to the linear 
order) from the requirement of the perturbation opera- 
tor being cancelled by this transformation. Therefore we 
transform the density operator according to 

nk = hk - [Ux,hk], (28) 

where Hi is the solution of the operator equation 

[Kx,H h ]=H* (29) 

under the constraint (1271) , that provides correct structure 
of the correct solution and ensures its uniqueness. 

Now solve Eq. (|2"9")h recalling the explicit form of the 
unperturbed Hamiltonian Eq (I10p and the perturbation 

operator Eq. (fl~7|) . We note that, obviously, B q _ q q _ ql = 
and B-q^qi qq i =0. Then symmetry arguments allow us 
to state that the coefficient B~ q ~ q q q does not change 
under any permutation of its arguments and B~ q , q = 

B~, q ~ q . After some algebra we obtain a set of equations 
for these non-zero coefficients, which in the matrix form 
reads as 



-3[(fc-g) 2 /4 4 
-3(<?74+l) 
-3(fc 2 /4+l) 




{k-qf 


-(P/4+1) 



k 2 


(k~q) 2 

-(g 2 /4 + i) 





k 2 

„2 



-[(fc-g) 2 /4+l] 



i B k-, 

k 

k-q q 



( "ID 

nmB q k _ q ~_ k 
\niuB h :i q q 



k q 



qk/2 
k{k-q)/2 
-q(k-q)/2 
~(k 2 - kq + q 2 )/8 



m 



The solution of the set of Eqs. (|30|) is 

Tlk-q -kq _ n lD 

3 ' 



B 



k-q q 



R q 

D k-q -k 



B k-q _ 



1 


'1 


i 




riiD 


4 ~ 


q(k- 


q). 


1 


'1 


1 




"ID 


4^ 


k{k- 




-5- 










a- 







(31) 
(32) 

(33) 
(34) 



Then we find the explicit form of the transformation 
Eq. (Ell) 



rik 



n k + 
1 



q(k - q) 



niT>4>k- q 4>q + 



Uk-qUq 



(35) 



III. TRANSITION MATRIX ELEMENT 



We transform the integrability-breaking interaction 
term Eq. flT3l by changing hk to hk, applying 
Eq. (|35[) . and retaining, according to the assumed or- 
der of approximation, only quartic terms. Then af- 
ter some algebra we obtain the matrix element (ko — 

h - k 2 , h, k 2 \Hni\ko), where \k ) = &[.J0), |fc - fa - 

k 2 , fa, k 2 ) = V ko _ kl _ k: p kl b\ 2 \Q) , and |0) is the vacuum of 
Bogoliubov elementary excitations. In a case of arbitrary 
initial numbers of the involved elementary excitations, 
the relevant matrix element 



(f|ff n i|in) = \J (nfci + l)(«fc 2 + 1 )(«fe -/c 1 -fc 2 + IK x 
(k - fa - k 2 , fa, k 2 \H n i\k ), (36) 



where 



in) = (n kl lnkJn ko -ki-k 2 } -nko } -) 1/2 x 
ft n k ft nk -*i-*2 £t " fe i £t "fc 2 m\ 

°fe u k -k l -k 2 °ki u k 2 l U /' 

|f> = [K+l)!K + l)!x 

(nk - kl -k 2 + l)!(n fco - l)!]- 1 / 2 x 



ft "fe - 

°k 



k — k 1 —k 2 fel k 2 



(37) 



|0>,(38) 



(|36)) enables us to evaluate the transition rates in a qua- 
sicondensate at non-zero temperature. 

After some algebra we obtain 



(fc - ki - k 2 , fei, k 2 \H ni \k Q ) = 

'-Mi 



12K 3 b . 



Lnn 



(39) 



is readily expressed through (fco — fci — k 2l fei, fel-ffnilfeo) 
and matrix elements of the bosonic field operators. Eq. 



where the dimcnsionless form of the matrix element is 



M k f° a ki k2kik2 = 5fe _fe 1 _fe 2 5fc 2 yk fcj + \J Sk -k 1 -k 2 Sk 1 yk„ k 2 + V Sk ± Sk 2 yk k a -k 1 -k 2 + 

\J Sk Sk 1 Zk -k 1 -k 2 k 2 + \J Sk Sk 2 Zko-k!-k 2 fei + \/ Sk Sko-k!-k 2 Zkt k 2 i (40) 

I 



z ---Jvm + (i-^)^- (42) 

The significance of our method to derive Eqs. (|39l 
can be understood at the best from comparison to 
a "naive" alternative way of derivation. Namely, we 
may express the density operator in Eq. (|5|) through 
the bosonic field creation and annihilation operators as 
h = tp^ip. Although there is no true condensate in the 
thermodynamic limit in ID, we may try, by integrating 
out slow variables [151 ] (up to the infrared cut-off mo- 
mentum fcipi), to express the bosonic field operator as 
i> = \/niT> + Sip, where Sip = £~ 1/2 J2\k\>k m ake tkz and 
the bare atomic operators a/j in the momentum represen- 
tation are related to the quasiparticle operators bk, b*_ k 
via the Bogoliubov transformation. The correct value of 
the rate of Beliaev and Landau damping in weakly inter- 
acting two-dimensional Bose gases (where true conden- 
sate is absent at any finite temperature) can be found 
in this way (23|, therefore it is not clear from the very 
beginning that this method fails for three-body colli- 
sions in ID. After renormalizing the two-body interaction 
strength and singling out the integrability-breaking terms 
terms one could obtain, to the lowest order in nij^a s an 
incorrect form for the integrability-breaking interactions 

H ni = : / dz (8ipi + SipfSipHip : (43) 

nm Jo 

(here symbol : : denotes normal ordering of the atomic 
field operators), that yields the result similar to Eqs. 
(SHI HO), but with y qq ,, Z qq , substituted by 



Z -' - 4v /V^ + 4 ' (45) 

respectively. As we see later, Eqs. (PHI l4"5t lead to a sin- 
gularity of the matrix element M k k °~ kl ~ k2 kl k2 at vanish- 
ing momentum transfer, thus signifying the inapplicabil- 
ity of the "naive" approach. 

To the end of this Section, we discuss the kinematics 
of the three-body process shown in Fig. [1] and the values 
of the matrix element (14U1) on the energy shell, where the 
energy conservation law requires 

£fe = efci + £fe 2 + e/co-fci-fcs ( 46 ) 

and the energy of an elementary excitation is given by 
Eq. (|14[) . Two momenta of the product elementary exci- 
tations (we denote them by fci and ko — k\ — k 2 , assuming 
for convenience |fci| < \ko — ki — k 2 \) have the same sign 
as fco, and the third excitation momentum has the op- 
posite sign, sgn(fc2fco) < 0. To be definite, we assume 
fco > 0. 

A. Phononic limit 

If feo ^ 1 then all relevant excitations are phonons 
having almost linear dispersion law. Only taking into 
account the cubic correction ss |fc| + |fc| 3 /8 helps us to 
find from Eq. (j46j) the relation between the momenta of 
the phonons: 

3 

k 2 feo(feo-fci)fei- (47) 

The matrix element (ES> calculated using Eqs. (|4Tl |42|) 
in the phononic limit is 
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FIG. 2: Feynman diagrams of the processes providing ther- 
malization of low-energy excitations (phonons). 



On the energy shell, where Eq. 
reduced to 



M 



ko—k\—k2 k\ k% 



1ZD holds, Eq. (HSJ is 



(49) 



V3k Q 



The matrix element (|49p is finite in the limit of vanish- 
ing transferred momentum, k± —¥ 0, unlike the incorrect 
"naive" estimation that relies on Eqs. (|4"4l |4"5")) and pre- 
dicts the matrix element diverging as l/k%. 

However, this process alone is suppressed for low- 
energy excitations (phonons) in trapped condensates. 
The reason is the infrared momentum cut-off imposed 
by a finite length £ of a trapped ID quasicondensate. 
From the condition \k\2 ^ 2ir/£ and Eq. (l47t we con- 
clude that this process lead to damping of phonons with 
k a £ 5£~ 1/3 (or, in usual units, fc > 5^ 2/3 £- 1 / 3 . In a 
87 Rb quasicondensate with the density nm = 50 /xm , 
length £ = 100 /mi, and the radial trapping frequency 
lj± = 2ir x 3 kHz the process shown in Fig. [1] is kine- 
matically allowed for fco > 3.5 x 10 4 cm -1 , which is very 
close to the crossover fc ~ 1/Ch between the phononic 
and particle-like parts of the Bogoliubov excitation spec- 
trum. 

The processes shown in Fig. [5J in contrast, have less 
severe kinematic restriction, ko > 2ir/£, and therefore can 
thermalize elementary excitations deeply in the phononic 
regime. 



B. Fast- particle limit 

In this limit ko ^> 1. As concerns the the product 
states, there are two distinct cases [T3|. All three product 
elementary excitations can be fast as well, in this case 
Eq.((46]) yields 



- I yjk* + 2fc fc 1 - 3kl (50) 




FIG. 3: Transition matrix element M = Ml 1 2 1 2 as a 

HQ 

function of the product state wavenumber fci. Units on axes 
are dimensionless. On the horizontal axis label we explicitly 
indicate scaling of fci to the inverse healing length. M is 
calculated on the energy shell, i.e. the energy conservation 
Eq. (|46|) holds. Solid lines: the main result of Sec. IIIII 
given by Eqs. (|40l - |42)| . Dashed lines: approximation by 
Eqs. (1441 1451) (singular at fci — > 0). Initial momentum fci = 
10 Ch 1 (thin lines) and 25 C^ 1 (thick lines). Fitting by Eq. 
(|49[) is practically indistinguishable from the solid lines and 
is displayed in the inset. 



The aforementioned ordering convention fc 2 < 0, < 
fci < ko — k\ — ki-, ka > now results in 



< fci < 2fc /3. 



(51) 



But there is another possibility, when two product el- 
ementary excitations are phonons, fci < 1 and 



fcp-1 
k + l 



fci. 



(52) 



Although the phase space corresponding to the latter case 
is relatively small, collisions with small transferred mo- 
mentum are bosonically amplified in a quantum gas with 
high thermal population of phononic modes, and their 
contribution to the fast particle damping can be signifi- 
cant fli ]. 

Eqs. (|40l - 142|) yield on the energy shell, i.e. for Eqs. 
(I5TJ1) and (|52"j) holding for fc x 1 and k 1 < 1, respectively 



K / ,k -k 1 -k 2 k 1 k 2 _ J 3- v /|fcifc 2 |/4, fci < 1 ,_„^ 
Mk « "I 3/2, fci»l ■ [b6) 

The formula interpolating between these two cases is 

3 



\yfk — fci— k 2 fci fc'2 

fco 



;V Sk ± Sk 2 



(54) 
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in full agreement with the used in Ref. [l4j assumption 
of correlation separation 

(ko — kx—kz ki k 2 \h 3 \k ) = 3(ki fc 2 |n 2 |0)(fco-fci-fc 2 |? : i|fco) 

(55) 

that singles out the fast particle contribution. 

In Fig. |3] we display the matrix element Eq. P0|) 
calculated in the case of fco 3> 1 with y q q > , Z q q i given by 
Eqs. (|4T1 142|) . and, for comparison, its value calculated 
from Eqs. (j43]-|45|). We see that the latter approach 
works only for asymptotically large k\, and the former 
one is in a very good agreement with the approximation 

(El- 
Approximation (|54[) breaks down out of the energy 
shell. This can be important for transient processes, with 
the energy uncertainty inversely proportional to the du- 



ration of the external driving of the system, like in the 
experiment [2J|. However, this problem (also related to 
the quantum Zeno and anti-Zeno effects [IE]) is to be 
considered separately from the kinetic equation, which 
follows from Fermi's golden rule and is considered in the 
next Section. 



IV. KINETIC EQUATION 
A. General approach 

Using Fermi's golden rule, we can readily write a ID 
kinetic equation that takes into account effects of bosonic 
amplification [cf. Eq. (|36l) ]: 



Of 



-7o 



dkidk 2 (5(e fc _ fcl -fc 2 + e kl + e k2 - e k ) 



Qi 



j^k — k 1 —k 2 fci fc s 



[(/fc-fci-fca + l)(/fci + l)(/fe 2 + 1)A- - fk-k 1 -k 2 fk 1 fk 2 (fk + 1)] + 
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dkxdk 2 S(e k + Cfcl + e k2 - tk+k 1+ k 2 ) M k k l\« 2 +k2 
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I 



(56) 



Here fk is the population of the elementary excitation 
mode with the ID momentum fc. The rate of the process 



is scaled to 70 = (72/7r)(K 3 b/n 



ID J 



or, in usual units, 



70 



18 2 / niDo; 
— wj_4 



(57) 



where l± = y^h/(muj±) is the radial trapping length 
scale. The first term in curly brackets in Eq. (|56p corre- 
sponds to processes shown in Fig. [1] with the wave number 
k assigned to the (single) incoming line in Fig.QJa) or the 
(single) outgoing line in Fig. Hfb). The second term de- 
scribes the situation of k being one of the three incoming 
quasiparticle momenta in Fig. [TJb) or one of the three 
outgoing quasiparticles in Fig. [Ha). The respective in- 
tegration ranges take into account the bosonic nature of 
elementary excitations and are defined in a way that pre- 
vents double counting of the same momenta of elemen- 
tary excitations. Qi 2 includes all fci, k 2 satisfying the 
conditions sgnfci = sgn(fc — k\ — k 2 ) — ~sgnk 2 = sgnfc, 
I fci| < \k — ki — k 2 \, < |fc2| < 00; Q' 12 consists of two 
regions: (i) sgn/c = sgnfci = — sgnfc 2 = sgn(A: — k\ — k 2 ), 
< \k\\ < 00, < \k 2 \ < 00, and sgn/c = —sgnfci = 
— sgnfc 2 = — sgn (fc— fci — fc 2 ), < |fci| < ki — k 2 \ < 00. 

In trapped quasicondensates, Eq. (|56[) applies in the 
local density approximation if discreteness of the ele- 
mentary exciation spectrum is not important, i.e., if 
\k 2 \ ^> 2ir/£. Practically, this means Eq. (|5r?|) applies 



in the fast particle damping regime for scattering events 
with \ki\ « |fc 2 | 2ir/£. In this case, competing pro- 
cesses shown in Fig. [5] should provide relatively small 
contribution to the thermalization rate, however, their 
detailed calculation is to be a subject of future work. 
A rough estimation (based on evaluation of phase space 
available for scattering of a fast particle with wave num- 
ber fc ^ 1 on other quasiparticles) predicts that the di- 
agrams from the first and second raws yield the rates 
smaller by a factor fc/niD and (fc/nio) 2 , respectively. 

It is easy to check that the stationary solution of Eq. 
(1561) is the Bose-Einstein equilibrium distribution 



1 



exp(e fc /T)-l 



(58) 



Here we measure temperature T in energy units (i.e., 
set Boltzmann's constant to 1). Eq. ((55)) does not con- 
serve the total number of elementary excitations, there- 
fore chemical potential in the Bose-Einstein distribution 
(l57l is zero. 

If the energy e k of the fast particle is much larger than 
both the temperature and mean interaction energy per 
particle then scattering with large transferred momen- 
tum populates initially empty particle-like modes, and 
the population of the fc-mode decreases exponentially 
with the decrement T dam P = \Z3ttj /4 [HG3 (note that 
pdamp (jpgg no t depend on fc). 
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It has been noticed (l4j that three-body collisions with 
small momentum transfer can give, due to bosonic am- 
plification, major contribution to the damping rate of a 
fast particle in a certain parameter range. In the present 
paper we make a step forward compared to the treat- 
ment of Ref. [lH by looking precisely to the dynamics of 
a bunch of fast particles in a ID quasicondensate. 

B. Fokker-Planck equation 

We consider a non-equilibrium distribution of elemen- 
tary excitations f k — f% + Sfk, where the perturbed part 
is small, dkSfk <C dkf^, and is localized in 
the momentum space around fco > 1 (the width of the 
fast particles bunch being Afc <C fco)- Such a distribu- 
tion can be created by means of Brag g sp ectroscopy [26j , 
perhaps, using higher-order processes [271 ] to obtain large 
values of ko. To be still in the ID regime and avoid scat- 
tering of atoms to radially excited states we must require 
the kinetic equation of the fast particles to be lees than 
2huj± (factor 2 appears here due to parity conservation 



[T2T [Hj]). A relatively small numbers of fast particles 
allows us to neglect heating of the lower modes in the 
course evolution. 

In what follows we take into account only collisions 
with small momentum transfer. Then we analyze the 
dynamics of 5f k on the time scale much shorter than 
1 /T damp (when we can neglect scattering with large mo- 
mentum transfer). A change of 6f k on such a short time 
seal will be a measure of the importance of the scattering 
events with the small momentum transfer, bosonically- 
amplified by f° ± ~ /| 2 » 1. 

The assumptions listed above enable us to linearize the 
kinetic equation ([551) and reduce it to the Fokker-Planck 
equation [28[. We neglect the momentum dependence 
of the kinetic coefficients (diffusion D and advection A 
within the narrow momentum around fco an d finally ob- 
tain 

l6f h =A^ + D^Sf kt (59) 



poo pQ 

A = 70 / dfci / dk 2 (fci + fc 2 ) 

JO J -00 

D = 70 / dkt / dk 2 (fci + k 2 f 



j^k -ki-k 2 ki k 2 



K / jk -k 1 -k 2 k 1 k 2 
ko 



<5(e fco ^fe 1 -fe 2 + £fc! + 6fe 2 - £fc )(/fc 1 + /fc 2 ), 



(60) 



S(e k0 ^ k2 + e kl + e k2 - e ko ) [/£ + (/ fc ° + (61) 
I 



Eqs. ([591 - [6"Tjl are written still in dimensionless vari- 
ables. In what follows we assume usual time and length 
units in Eqs. ([59p and evaluate A and D in two limiting 
cases. First, for T < giD^iD we obtain 



14.6 sgn ko 



70 



T 



D 



2.45 



70 



fc o| 2 C/j VSlD^lD 

T 



(62) 
(63) 



where 70 is given by Eq. ([57[) . 

If, on the contrary, gmniD <T^ 2huj± then bosonic 
amplification of scattering to free-particle-like mode be- 
comes significant and we obtain 



A ss 14.7sgnfc 



On 



\ko\ 2 Q \giun 



T 



D 



10.8- 



7o 



T 



k o\ 3 Cl V^ID^ID 



'ID 



5/2 



3/2 



(64) 
(65) 



The decrease of the kinetic coefficients A and D with 
growing fco is explained by two observation. First, the to- 
tal momentum transfer decreases as fco increases but the 
low-energy excitation momentum k\ is kept fixed (and 



k 2 is related to fco a- n d fci via energy and momentum 
conservation): fci + fc 2 ~ |fci|(4CiT 2 + fc 2 ) 1 ^ 2 /fco- Second, 
integrating in Eqs. (pOl l6"Tj) the delta-function of the en- 
ergy difference over fc 2 yields another factor 1/fco- 

Assume that the initial distribution of fast particles 
is Gaussian, centered at fco and with r.m.s (determined, 
for example, by the duration of the Bragg pulse [26|) 
Afc. The well-known solution [H| of Eq. $5§§ with the 
Gaussian initial condition is 5f k (t) cx exp[— (fc — fco + 
At)7(2Afc 2 + 4L>t)]/(2Afc 2 + 4Dt)^ 2 . We define the 
typical time scales tA and tjj as times when shift of the 
maximum of Sf k or increase of its width, respectively, 
become comparable to the initial width Afc: 



t A = Ak/A, t D = Afc 2 /(2L>). 



(66) 



Bosonically-amplified three-body collisions with small 
momentum transfer are practically important, if they 
modify Sf k faster, than scattering with large momentum 
transfer to initially empty modes, i.e., if at least one of 
two inequalities 



T damp t A 



<1, 



p 



1P *D < 1 



(67) 



is satisfied. Note, that the number of scattering events 
per unit time per one fast atom can be significantly larger 
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40 60 80 100 

T(nK) 

FIG. 4: Typical times of shift (£a, solid lines), diffusion (tD, 
dashed lines), and scattering with large momentum transfer 
(t F = i/r damp , dot-dashed lines) for 87 Rb quasicondensate 
with niD = 40 /an" 1 (thin black lines), 60 /im -1 (medium 
gray lines), and 80 /xm -1 (thick black lines) as functions of 
temperature. cj_l = 2-7T x 3 kHz. Initial distribution of fast 
atoms is Gaussian with the mean ko = \l 2muj±/h = 0.7 x 
10 5 cm -1 and the r.m.s. Afc = jko- 

than r damp sa 1.36 70, but their influence on the dynam- 
ics of the distribution 5fk may be still small. 

Fig. g] illustrates t A , t D juxtaposed with t F = l/^ dam P 
for experimentally realistic system parameters. One can 
see that, at experimentally feasible densities, first the 
momentum-shift time t A and then the diffusion time to 
become shorter than the typical time t F , as temperature 
increases. However, we should always keep in mind the 
limitation T < hui± necessary for the ID regime. 

Now we can compare our criterion for the significance 
of the small-transferred- momentum scattering events and 
that proposed by Tan, Pustilnik and Glazman [l4[ ■ In the 
latter case, the small-transferred-momentum scattering 
events were considered important if the total collision 
rate per atom far exceeded the value of r damp . This 
corresponds, in the limit 

fcoCh > 1, (68) 

to temperatures T > giDnmikoCh) 1 ^ 2 ■ We show, that 
this condition is not sufficient for experimental observa- 
tion of the effect of these collisions, and require Eq. (foT)) 
instead. The shift of the mean momentum of fast parti- 
cles becomes significant at 

?i 2 fc 4/3 Afc 2 / 3 

T > 0iDn 1D (fcoa) 4/3 (AfcOO 2/3 = ° , (69) 

m 

and large diffusive spreading of the momentum distribu- 
tion of fast particles requires 

T > 9i D n 1D (fcoO0 6/5 (Afca) 4/5 = ^ fc ° 6/5Afc4/5 - (70) 



Note that Eq. (|68|) does not hold for standard two- 
photon Bragg spectroscopy (26l |. which can provide 
koCh ~ 5. In this case, a rough estimation yields, instead 
of Eq.flMlEOl), T £ SidUd (cf. Fig. H. 



V. DISCUSSION AND CONCLUSIONS 

Eq. (l67|) is not the only restriction that applies in 
a real experiment. In practice, we always need to take 
into account the finiteness of the system. The longitu- 
dinal trapping determines the finite size of the ultracold 
atomic cloud [29j]. This size can be associated, by or- 
der of magnitude, with the parameter t introduced in 
Sec. IIII Al The longitudinal trapping potential itself 
lifts, in principle, the integrability of the Lieb-Liniger 
model, however, we expect this non-integrability effect to 
be small [3(| ■ More important is anharmonicity (of the 
longitudinal trap itself or induced by the atomic mean- 
field interactions) that induces dephasing and thus lim- 
its the time of observation of the integrable dynamics 
to few oscillation periods Q- The most desirable case 
takes place if the bosonically-amplified three-body elas- 
tic processes manifest themselves after a single passage 
of fast atoms through the cloud, i.e. on the time scale 
t = m£/(hko). The Bragg pulse that produces the fast 
atoms should be also much shorter than To. In contrast 
to the Fourier-limited width of the excitation spectrum in 
the frequency domain, the momentum distribution width 
Afc rj mT/(2h 2 n 1D ) of the fast particles is determined 
by thermal fluctuations of the phase of the quasicon- 
densate [3l|. Regarding a longitudinally trapped qua- 
sicondensate, we will use the notation nxo f° r its average 
linear density. Although the momentum distribution is 
Lorentzian (3j |32| , we still use Eq. (|6"6"|) , which is derived 
in the case of a Gaussian, for rough estimations. We as- 
sume that, after the Bragg pulse, fast atoms are left to 
interact with the rest of the quasicondensate for the time 
^int ~ t~q. Then the trap is suddenly switched off, and 
the atomic cloud expands ballistically. The atomic posi- 
tion distribution in the asymptotic regime of long time 
of flight corresponds to the in-trap momentum distribu- 
tion [2g. We assume that the influence of the bosoni- 
cally amplified three-body collisions with small momen- 
tum transfer is unambiguously detectable if the relative 
shift of the momentum of fast stoms due to three-body 
elastic scattering is of about 50%, i.e. At int s» 0.5 fco- To 
be definite, consider an atomic cloud of size I w 100 //m. 
Initial wave number fco = 0.7 x 10 5 cm -1 of the fast par- 
ticles corresponds to To ~ 20 ms. Taking i; nt = 10 ms, 
we find that the bosonically-amplified three-body scat- 
tering is detectable at T > 110 nK for niD = 80 /zm _1 
and T > 140 nK for ?iid = 60 /im _1 . The linear density 
niD = 40 /im^ 1 is too small, and the slowing down of 
fast particles can not be seen in the temperature range 
corresponding to the ID regime. Note for all these den- 
sity values r damp tj nt < 1, i.e. three-body scattering into 
empty modes (with lage momentum transfer) is too weak 



10 



to be detected in a singe passage of fast atoms through 
the quasicondensate. 

An alternative (to Bragg spectroscopy) way of produc- 
ing fast atoms in twin-atom beams by external driving 
of the radial excitations of ultracold atoms in an elon- 
gated trap has been recently demonstrated by Biicker et 
al. 

Finally, we discuss the relation between the theory 
of ultracold atomic systems in the quasi- ID regime and 
the quantum Newton's cradle experiment Q that has 
brought a new attention to the problem of the integra- 
bility breakdown. Thcrmalization has not been observed 
in that experiment, and only lower limits to the ther- 
malization time have been obtained for different (strong 
and intermediate) interaction strengths [7(. This fact is 
in a quantitative agreement with the estimations of the 
thermalization rate [l4[ (see also Ref. [l3j]) taking into 
account the effect of the atomic correlations [the zero- 
distance two-body correlation function 52(0) < 1] in the 
case of strong repulsive interactions. The damping time 
1/ [7032(0)] occurs to be too long to be measured. Also a 
small size of the colliding atomic clouds puts a relatively 
high infrared momentum cut-off (2ir/£ > 0.1 /1m -1 ) that 
precludes scattering events with small momentum trans- 
fer. As we mentioned before, non-integr ability caused by 
the longitudinal harmonic confinement is also too weak 



to be detected [30J . Another possible mechanism of the 
non-integrability, associated with the tunnel coupling be- 
tween adjacent waveguides in a 2D optical lattice, can be 
relevant only for optical lattices much weaker than that 
applied in Ref. [7| . To observe thermalization of bosonic 
atoms in a quasi-lD waveguide and, in particular, to test 
the theory developed in the present paper, one has to 
work at larger atom numbers that provide both higher 
atomic linear densities and longer system sizes than in 
the quantum Newton's cradle experiment Q- 

To summarize, in this paper we analyzed effective 
three-body collisions (mediated by virtual excitations of 
the radial degrees of freedom) in a ID quasicondensate of 
ultracold bosonic atoms. We calculated the matrix ele- 
ment for the decay of a single elementary excitations into 
three elementary excitations. We stress that the obtained 
expression is non-divergent at small momentum transfer. 
We derived a kinetic equation governing the damping of 
fast particles in a quasicondensate and its Fokkcr-Planck 
limit that accounts for scattering into thermally popu- 
lated modes with small momenta. We demonstrate that 
the latter process can be observed experimentally. 

The author thanks J. Armijo, I. Bouchoule, L. I. Glaz- 
man, H. Grosse, S. Manz, J. Schmiedmayer, and J. Yng- 
vason for helpful discussions. This work is supported by 
the FWF (project P22590-N16). 



[1] E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963); 
E. H. Lieb, Phys. Rev. 130, 1616 (1963). 

[2] H. B. Thacker, Rev. Mod. Phys. 53, 253 (1981). 

[3] T. Kinoshita, T. R. Wenger and D. S. Weiss, Science 305, 
1125 (2004); O. Morsch and M. Oberthaler, Rev. Mod. 
Phys. 78, 179 (2006). 

[4] R. Folman, P. Kriiger, D. Cassettari, B. Hessmo, T. 
Maier, and J. Schmiedmayer, Phys. Rev. Lett. 84, 4749 
(2000); J. Fortagh and C. Zimmermann, Rev. Mod. Phys. 
79, 235 (2007). 

[5] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

[6] Note that it is practically sufficient to have temperature 
and the mean-field interaction energy slightly below (or 
even comparable to) Juj± to ensure the ID regime, as 
demonstrated in many experiments: S. Manz, R. Biicker, 
T. Betz, Ch. Roller, S. Hofferberth, I. E. Mazets, A. 
Imambekov, E. Demler, A. Perrin, J. Schmiedmayer, 
and T. Schumm, Phys. Rev. A 81, 031610 (2010); P. 
Kriiger, S. Hofferberth, I. E. Mazets, I. Lesanovsky, and 
J. Schmiedmayer, Phys. Rev. Lett. 105, 265302 (2010); 
T. Betz, S. Manz, R. Biicker, T. Berrada, Ch. Roller, 
G. Razakov, I. E. Mazets, H.-P. Stimming, A. Perrin, T. 
Schumm, and J. Schmiedmayer, Phys. Rev. Lett. 106, 
020407 (2011). 

[7] T. Kinoshita, T. Wenger, and D. S. Weiss, Nature 440, 
900 (2006). 

[8] S. Hofferberth, I. Lesanovsky, T. Schumm, A. Imam- 
bekov, V. Gritsev, E. Demler, and J. Schmiedmayer, Na- 
ture Physics 4, 489 (2008). 

[9] L. Salasnich, A. Parola, and L. Reatto, Phys. Rev. A 65, 
043614 (2002). 



[10] A. Muryshev, G. V. Shlyapnikov, W. Ertmer, R. Seng- 
stock, and M. Lewenstein, Phys. Rev. Lett. 89, 110401 
(2002). 

[11] L. Rhaykovich and B. A. Malomed, Phys. Rev. A 74, 

023607 (2006); L. Salasnich and B. A. Malomed, Phys. 

Rev. A 74, 053610 (2006). 
[12] I. E. Mazets, T. Schumm, and J. Schmiedmayer, Phys. 

Rev. Lett. 100, 210403 (2008). 
[13] I. E. Mazets and J. Schmiedmayer, New J. Phys. 12, 

055023 (2010). 

[14] S. Tan, M. Pustilnik, and L. I. Glazman, Phys. Rev. Lett. 
105, 090404 (2010). 

[15] V. N. Popov, Functional Integrals in Quantum Field The- 
ory and Statistical Physics (Reidel, Dordrecht, 1983). 

[16] C. Mora and Y. Castin, Phys. Rev. A 67, 053615 (2003). 

[17] M. Lewenstein and L. You, Phys. Rev. Lett. 77, 3489 
(1996). 

[18] S. T. Beliaev, Sov. Phys. JETP 34, 299 (1958). 
[19] W. V. Liu, Phys. Rev. Lett. 79, 4056 (1997). 
[20] P. C. Hohenberg and P. C. Martin, Ann. Phys. (N. Y.) 
34, 291 (1965). 

[21] A. Ishihara, Statistical Physics (Academic Press, New 
York, 1971). 

[22] J. J. Sakurai, Modern Quantum Mechanics (Addison- 

Wesley, Reading MA, 1994). 
[23] M.-C. Chung and A. B. Bhattacherjee, New J. Phys. 11, 

123012 (2009). 

[24] N. Bar-Gill, E. E. Rowen, G. Rurizki, and N. Davidson, 

Phys. Rev. Lett. 102, 110401 (2009). 
[25] A. G. Rofman and G. Rurizki, Nature 405, 546 (2000). 
[26] D. M. Stamper-Rurn, A. P. Chikkatur, A. Gorlitz, S. In- 



11 



ouye, S. Gupta, D. E. Pritchard, and W. Ketterle, Phys. 
Rev. Lett. 83, 2876 (1999); J. Steinhauer, R. Ozeri, N. 
Katz, and N. Davidson, Phys. Rev. Lett. 88, 120407 
(2002). 

[27] E. E. Rowen, N. Bar-Gill, and N. Davidson, Phys. Rev. 

Lett. 101, 010404 (2008) 
[28] C. W. Gardiner, Stochastic Methods (Springer, Berlin, 

2009). 

[29] J. N. Fuchs, X. Leyronas, and R. Combescot, Phys. Rev. 
A 68, 043610 (2003); F. Gerbier, Europhys. Lett. 66, 771 
(2004). 

[30] I. E. Mazets, arXiv: 1011.0713 (accepted for publication 
in Eur. Phys. J. D). 



[31] S. Richard, F. Gerbier, J. H. Thywissen, M. Hugbart, 
P. Bouyer, and A. Aspect, Phys. Rev. Lett. 91, 010405 
(2003). 

[32] The shape and the width of the momentum distribution 
of atoms decoupled by the Bragg pulse [3l|] are in a good 
agreement with their estimation from the Fourier trans- 
form of the first-order spatial correlation function of a 
weakly-interacting quasicondensate [TH. [T^]. 

[33] R. Biicker, J. Grond, S. Manz, T. Berrada, T. Betz, Ch. 
Koller, U. Hohenester, T. Schumm, A. Perrin, and J. 
Schmiedmayer, arXiv: 1012.2348. 



